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Abstract 


The design of an experiment to measure skin friction and 
turbulent boundary-layer characteristics at Reynolds numbers 
exceeding 1 x 10 9 is described* The experiment will be conducted 
in a zero-pressure-gradient flow on a flat plate in the National 
Transonic Facility (NTF). The development of computational codes 
to analyze the aerodynamic loads and the blockage is documented* 
Novel instrumentation techniques and models, designed to operate 
in cryogenic environments, are presented. Special problems 
associated with aerodynamic loads, surface finish, and hot-wire 
anemometers are discussed. 
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List of Symbols 


Roman 

C f Local skin friction 
C F Average skin friction 
C u Lift coefficient at M=0 

C L Lift coefficient at operating Mach number 
C p Pressure coefficient 
c Chord, 5.2 m 
H Shape factor, S*/B 
k Allowable roughness height 
L Lift 
M Mach number 
m /1-M 2 

q Dynamic pressure 
Re Reynolds number based on x 
Re / Unit Reynolds number [m‘ l ] 

S Area of plate, 12.5 m 2 
T Temperature 
U Freestream velocity 
v* Shear velocity, /tan,/p 

Greek 

a Angle of attack (see Fig. 2) 

A Flap angle (see Fig. 2) 

A C p Pressure difference across plate 
Displacement thickness 
B Momentum thickness 
v Kinematic viscosity 
P Density 

t * Wall shear stress, C f pU 2 /2 
Subscripts 

0 Conditions at M=0 

1 Extreme test condition 
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I. Introduction 

Knowledge of the turbulent boundary-layer characteristics of 
high-speed flows is important for the design of any modern air- 
transport vehicle. The increase in size of these vehicles has 
pushed the aerodynamic design requirements to regions of higher 
chord Reynolds numbers, Re c , that approach 1 x 10’. For example, 
Re c = 60 x 10* on the 747 wings, Re c = 5 x 10® on the C5 fuselage, 
and Re c = 1 x 10 ? on the B-l fuselage. Submarines routinely have 
Reynolds numbers that exceed 1 x 10’. At the present time, these 
Reynolds numbers can be a factor of 6-100 times greater than 
those obtainable by existing facilities. 1 Moreover, there is a 
dearth of fundamental data on skin friction and velocity profiles 
at these Reynolds numbers for even the case of a zero-pressure- 
gradient, flat-plate flow. 

Figure 1 shows the effect of Reynolds number on flat-plate 
skin friction in the high-Reynolds-number range. For 

convenience, only those data of Winter et al. 2 are included. 
Other data at lower Reynolds numbers are 25-40 years old. The 
empirical curves 3 ' 6 show divergence when extended to the Reynolds 
number range of interest but not surprisingly since they were 
based on lower Reynolds number data (see White 7 for additional 
empirical models). This situation indicates further 

complications when estimates are made for the case of roughness 
or non-zero pressure gradient. 

With the anticipated operational use of the National 

Transonic Facility (NTF), we have the capability to do research 
and development in high-Reynolds-number flows. Of particular 
interest are data essential for evaluating turbulence theories 
and estimating high-Reynolds-number drag. Besides basic skin- 
friction data, many other unknowns exist. For example, the 

nature of large-scale structures within the turbulent boundary 
layer at high Reynolds numbers is uncertain. Indeed, their very 
existence in the Reynolds-number range of 10 ? can be questioned. 
Ultimately, the role of wall roughness, forward and backward 

facing steps, cavities, protuberances, and pressure gradients 
must be evaluated. Moreover, extensions to high Reynolds numbers 
of the basic results of turbulent-drag-reduction programs must be 
carried out. 

The zero-pressure-gradient, flat-plate flow is the intrinsic 
part of these more complicated flows. Therefore, a careful flat- 
plate experiment at high Reynolds numbers serves to establish a 
bench mark for future experimental work and theoretical modeling. 
This bench mark will extend the data base upon which future work 

will be done. The initial work will concentrate on basic skin- 

friction data with additional hot-wire measurements to explore 
the boundary-layer velocity profiles. 
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II. Test Conditions 


The general description and capabilities of the NTF are 
given by Fuller. 8 Therefore, only a few basic features are 
presented here. The NTF is a cryogenic, continuous-flow, high- 
pressure wind tunnel capable of operating at unit Reynolds 
numbers of up to 450 x 10 6 nr 1 . The test section is 2.5 m square 
and 7.6 m long with a slotted-wall configuration. Moreover, the 
initial configuration will have all of the slots closed. The 
test-section geometry can also be varied by slightly deflecting 
the upper and lower walls about a pivot point at the test-section 
entrance. This will be used to balance the displacement- 
thickness effects on the upper and lower walls during the zero- 
pressure-gradient tests. 


The temperature range of the facility extends from 100*K to 
340*K. In a range of Mach numbers around 0.7, chord Reynolds 
numbers can exceed 1.5 x 10’ on a model 4 m long. The high- 
Reynolds-number , high-dynamic-pressure test conditions occur for 
the low-temperature range of the facility and present a special 
set of problems with regard to models and instrumentation. 
Therefore, the test design uses the following conditions as 
extreme goals. 


Extreme Test Conditions 


Maximum 

Maximum 

Minimum 

Maximum 

Maximum 

Maximum 


Mach number Mj = 0.7 

dynamic pressure: qj = 2atm = 2 x 10 5 Pa (4300 lb/ft 2 ) 

temperature = Tj = 100'K (-270'F) 

unit Reynolds number: Re/ = 450 x 10 £ /m (137 x 10 £ /ft) 

wetted length: 4.1 m 

Reynolds number: 1.85 x 10 9 


The corners of the test matrix can be defined using Fuller 8 . 
For non-cryogenic operation, say air at 250*K, it is possible to 
achieve Re = 0.6 x 10 9 at M = 0.7. For 1 ow— Mach— numb er 

operation, say M = 0.3, one obtains Re = 0.7 x 10* at T = 100'K. 
At M = 0.7 and T = 100'K, Re = 1.85 x 10 9 can be achieved without 
having locally supersonic flow. Thus one can see that 

measurements can be made independent of either Mach number or 
static temperature. Thus, the role of Mach number and static 
temperature (although not considered to be important) will be 
evaluated along with Reynol ds-number effects. 

It has been established 9,10 that it is feasible to conduct 
experiments of this kind in a cryogenic environment. Kilgore et 
al . 9 did low Mach-number experiments on laminar boundary layers 
where the temperature was varied while the Reynolds number was 
held constant and observed no important differences. Adcock and 
Johnson 10 did computational simulations of turbulent boundary 
layers and showed that C f was within ±0.5% when real-gas effects 
were included at cryogenic temperatures and pressures 
characteristic of NTF. Of course, one of the tenets of this 
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program is to establish the validity of such computational models 
in these environments. 


Ill- Model 

The proposed experiment consists of a vertically-mounted 
flat plate that spans the tunnel as a splitter plate. In this 
configuration, the flat-plate flow is contaminated by turbulence 
and spanwise variations generated at the leading-edge/floor 
junction. The spreading angle is in the range of 11* to 15*, 
therefore, one can obtain about 4.6 m of uncontaminated flow. 
The pressure recovery from the trailing-edge flap moves upstream 
and begins to affect the flow at x = 4.1 m for a 5.2 m plate. 
Therefore, we make the conservative estimate that the maximum 
useful chord will extend to x = 4.1 m. 

The flat plate is 50 mm (2") thick with an overall 5.2 m 
chord (including leading edge and flap) and 2.5 m span. It is 
made of an aluminum skin, aluminum honeycomb construction. All 
parts and fixtures of the entire test assembly are made of 6061- 
T6 aluminum in order to avoid problems of differential thermal 
expansion. The skin thickness is 4.8 mm (3/16") consistent with 
strength requirements. The honeycomb is perforated and is 41.3 
mm (1-5/8") deep with a 6.4 mm (1/4") cell size. An epoxy 
certified to operate at 100*K is used. A number of special 
requirements have been developed to assure ourselves of achieving 
the test goals, and these are described in the subsequent 
sections . 


A. Surface Waviness 

Special care must be taken to control the surface waviness. 
The allowable amplitude A, of departures from flatness over a 
distance D, is given by A/D < 0.002. This criterion restricts 
the error in the skin-friction measurement to less than \% and 
follows from the work of Lin et al. u Because of this 
restriction, a large autoclave (3m x 6m) is used in the assembly 
of the model. 


B. Surface Roughness 

Because of the thin boundary layers, a very smooth and 
polished surface is required. A very conservative estimate 12 of 
the maximum allowable roughness height, k, is given by 

k = /50/C f /Re/ (1) 

where C f is the skin-friction coefficient and Re / is the unit 
Reynolds number. Equation (1) is derived by assuming that a 
surface is hydraulically smooth if all of the protuberances are 
contained within the laminar sublayer. This is true for v*k/v<5 



Final Report: Design of High-Reynolds 

Number Flat-Plate Experiments in the NTF 


Page 6 


where v* is the shear velocity /tan „/p. For the extreme test 
conditions, k =17 microinches over the test portion of the 
plate. Therefore, the plate surface-roughness criterion is 14 
microinch absolute and 12 microinch rms. Thus, careful hand 
polishing of the surface is required. 

C. Leading Edge and Trailing Edge 

The leading edge is a symmetric ellipse with a 12:1 ratio of 
major to minor axes. Its chord is 0.3 m and is made from a solid 
piece of aluminum. This configuration gives a good stagnation 
point and a tangent to the flat plate with only a small pressure 
gradient in the leading-edge region. It will be shaped on a 
numerically-controlled mill to hold the spanwise nonuniformities 
near the leading edge to the ±0.001 inch, after which it will be 
polished. A small strip on the leading edge will be left 
unpolished, and this should be sufficient to fix the transition 
location. 

A trailing-edge flap is used to control the position of the 
attachment line on the leading edge by balancing the blockage 
between the two sides of the p,late. The position of the flap is 
mechanically actuated to provide the desired flow conditions at 
different freestream velocities. 

Figure 2 shows the basic plate/flap configuration inside the 
tunnel. The plate is off-center and over the third slot from the 
far wall. The plate is at a slight negative angle of attack 
because of displacement-thickness blockage. The flap is adjusted 
to minimize the overall circulation. 


IV. Flow-Field Calculations 

A number of tasks are defined under flow-field calculations 
that must be carried out before model fabrication and testing 
begin. They are: (A) allowable pressure gradient for zero- 
pressure-gradient flow, (B) wall interference and blockage, and 
(C) aerodynamic loads on the model. 

The basic tools for the calculations are the panel-method 
code of Bristow 13 and the turbulent boundary-layer code of 
Anderson and Lewis. 14 These are used separately and in a 
combined version. 15 


A. Allowable Variation in C p 

In order to establish a criterion for zero-pressure-gradient 
flow, we assume that a local pressure gradient must not change 
the local skin friction by 1 % over the zero-pressure-gradient 
case. The momentum-integral equation 12 is usually written in the 
following form 
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xjp = (25 + «*)$*- + U 2 ^ (2) 

where is the wall shear stress, U(x) is the freestream 

velocity, 5 is the momentum thickness, and 5* is the displacement 
thickness. We need a relationship between the skin-friction 
coefficient and the gradient of the pressure coefficient so by 
using the definitions of C f , C p , and the shape factor H, the 
momentum-integral equation (2) can be rewritten as 

Cf = ~(2 + 3 (1 - C (3) 


For a zero pressure gradient 


C? = 


oAil 

dx 


(4) 


where superscript 0 denotes the zero-pressure-gradient reference 
point. By introducing a pressure gradient dC^ /dx into Equation 
(3), the other parameters change as C f = C| + CJ , 5 = 5° + O ' , 

H = 11°+ H / . Then the change in skin friction for a given 
gradient in C p can be written as 


AC, 


C2 


[-(2 + H«) 


dC^ 1_ d5» . 1 d5« s 

dx 50 dx ' '5 ( ' dx ' 


(5) 


where AC f = - C$. Equation (5) relates the change in skin 

friction as a function of the pressure gradient and the gradient 
of the perturbation of the momentum thickness. 


It is necessary to develop means for evaluating d5//dx. 
This is done using an approximation based on Buri’s method 12 for 
boundary layers with pressure gradients (see pp. 678-679 of 
Schlichting 12 ) . It can be shown after considerable algebra that 
to within the order of ( 0 . 0 1 ) Re^*^ 4 : 


d5/ 

dx 


dx 


Where A is a positive constant* of order one. Then Equation (5) 
is rewritten as 


AC f ,2 + H* , d C/ 1 d 5° . 

C° L 2 J dx ' V dx } 


( 6 ) 


In terms of an allowable 
expressed as 

dCC . 

dx ' allowable 


AC f /CJ , the allowable dC^/dx can be 


C|_ 1 , ac m 

5° 2+H°-2A ^ C f ° ; 


(7) 


*In terms of Buri’s method, A=2n/5 where n=3.94 and is defined in 
Equation 22.15 of Schlichting. 12 
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The constant, A, can be determined using selected boundary-layer 
calculations using Ref. (14) with Buri’s approximation. Thus 
Equation (7) (whose final form was suggested by Nagib I£ ) can be 
used as an approximate tool for estimating the allowable pressure 
gradient for an allowable error in C f . Moreover, Equation (7) 
defines a means for simplifying the number of turbulent boundary 
layer calculations by permitting a local evaluation with zero 
pressure gradient results. For the Reynolds number range of 10 9 , 
calculations show that 


dx ' allowable 1 M C? ' 


( 8 ) 


For the conditions of the present experiment, given an allowable 
error in C f of ±1%, a conservative allowable pressure gradient is 
dC^/dx < ± 2 x lO'Vcm at 3/4 chord. 


B. Wall Interference and Blockage 

A large number of turbulent boundary-layer calculations were 
made under different flow conditions within the test matrix. 
Based on the criterion given in Equation (8), it appears as 
though a plate orientation of approximately -0.1* angle of attack 
will account for displacement-thickness effects over the model. 
(See Figure 2 for plate orientation.) However, in order to 
establish this more accurately and to assure ourselves that wall 
liners and wall modifications are not required, the inviscid 
code 13 and the viscous code 14 are combined into a single blockage- 
calculation code. 15 

The C p distribution of the tunnel (beginning in the settling 
chamber) and of the plate (for some initial orientation of Figure 
2) are calculated first. These results are used to calculate the 
boundary layers on the sidewalls and on the plate. The shapes of 
the wall and of the plate are changed by adding the displacement 
thickness and the inviscid code is rerun. Only the blockage of 
the sidewalls is used since the positions of the floor and 
ceiling surfaces are adjusted independently for zero pressure 
gradient. This process is repeated until convergence. The angle 
of attack of the plate is varied in order to flatten the C p 
distribution to almost zero dp/dx, and the flap is adjusted to 
prevent separation on the test side of the leading edge. The 
inviscid calculations are done for M=0 and the C p distribution is 
scaled using the Prandt 1-G lauret rule where 

C, = C p0 ( 1-M 2 ) ' 1/2 =C |)0 /m (9) 

The boundary-layer code is run at M=0.7 considering this to be 
the most extreme condition. There are no locally supersonic 
regions in the flow at this Mach number. 

Figure 3 is a plot of -C Dfl vs x/c on the test side of the 
plate for an angle of attack of a=-0.08 and a flap angle of A=- 
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1.0. Between x/c = 0.15 (x=.78) and x/c=0.8 (x=4.1m), dC^/dx < 
1.0 x 10'Vcm which meets the requirements given in the 

discussion after Equation (8). The following table gives the 
values of the pressure gradient at M=0.7 for different chord 
locations . 


x/c 

dC p /dx x 10« 

0.117 

3.71 

0.161 

1.13 

0.214 

0.71 

0.390 

0.13 

0.507 

-0.001 

0.683 

-0.04 

0.742 

-0.36 

0.800 

-1.06 

0.860 

- 5.00 


The normalized numbers in the second column can be interpreted as 
dC p in X per meter. 

Thus it is shown that a zero-pressure-gradient flow can be 
established with a uniform angle of attack of the plate without 
resorting to sidewall modifications. This is known to be true 
for incompressible flows and is possible here provided the Mach 
number does not exceed 0.7. These results will be confirmed by 
conducting a scale model test in the Diffuser Flow Apparatus 17 at 
NASA Langley Research Center. The slightly negative C p0 shown in 
Figure is accounted for in the boundary-layer measurements by 
determining the virtual leading edge. 


C. Aerodynamic Loads 


Figure 2 shows the model configuration 
inside the tunnel. The conventions for the 
angle of attack, and flap are also shown, 
drawing. The average lift forces scale with 
rule and are given by 


as it is placed 
directions of lift, 
This is not a scale 
the Prandtl-Glauret 


L = C L( ,qj S ( 1-Mf )‘ I/2 (10) 

where C u is the lift coefficient determined by integrating, 
along x/c, the difference of the pressure coefficient, AC p0 , on 
each side of the plate. 


Figure 4 gives AC p0 along the model for a = -0.08* and A=- 
1.0*. The angle of attack produces a zero dp/dx over the test 
side of the plate and a pressure drop along the back side of the 
plate. Under these conditions, C u = #0.030 (the lift is 
positive and in the direction shown in Figure 2), and the total 
lift under extreme conditions is 1.1 x 10 5 N (25,000 lbs.) which 
represents an average loading of 8.5 kPa (1.2 psi). The loads in 
the leading edge and the flap are each 11,000 N (2500 lb.) 
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assuming an average C Lc of 0.05. The presence of the tunnel 
walls essentially doubles the lift over the no-wall case (an 
external flow). For example, with zero flap deflection we find 
that C u % 4 no.. 

The flap position can strongly effect the overall lift on 
the plate. For example, A=-1.5* gives C u =+0.054 while A=-0.5* 
gives C u =-0.0072. By setting the flap angle with slightly 
positive C u , the position of the attachment line on the leading 
edge does not separate the flow on the test side and compensates 
for extra blockage on the non-test side of the model due to 
instrumentation. The flap position is adjustable during the test 
because there is no way to account for this blockage. All load 
calculations are done at M=0.7 although C p0 and C Lo are used for 
reference. 

A computer code is written to calculate the deflections of 
the model using plate theory under load conditions typical of 
Figure 4. It assumes that the web of the honeycomb offers no 
strength and that all of the load is taken in the skin. The 
plate loading of 1.2 psi is well within the strength 
characteristics of the model. What is necessary, however, is to 
restrict the plate deflections so as not to violate the surface 
waviness criterion of section III. A or the allowable pressure 
gradient criterion of Equation (8) in section IV. B. It is 
assumed the plate has six supporting struts on the non-test side 
located at 1/4, 1/2, 3/4 chord and 1/3, 2/3 span. Small wide- 
flanged beams extend along the chord at these span locations. 
Fixed-end supports at the floor and ceiling are also assumed. 
Maximum deflection occurs mid-span at x/c=0.63 where AC p0 =O.O3. 
For the model configuration of section III, the maximum deflection 
is found to be 0.3 mm (0.012"). This is well with the surface 
waviness criterion. 


V. Instrumentation 

The model is instrumented with static-pressure ports, shear- 
stress gauges, thermocouples, strain gauges, and an 
accelerometer. Single location boundary-layer measurements are 
conducted with a total-pressure tube and a hot wire mounted on a 
variable jack set through the shear-stress-gauge housing. 

A. Pressure and Temperature Measurements 

There are a total of 107 static-pressure ports distributed 
along the plate as shown in Figure 5. The pressure ports are 6.3 
mm diameter inserts with a 0.3 mm orifice and are inserted into 
the plate before the surface is finished. The locations of the 
fifteen thermocouples are also shown in Figure 5. These are also 
mounted on inserts that fit into the plate. With these arrays, 
the wall conditions can be accurately monitored. Although the C p 
and A C p are not expected to exceed 5%, the high dynamic pressures 
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produce maximum 
order to satisfy 
need to resolve 
test conditions 
These values are 


differential pressures of around 0.1 atm. In 
Equation (8) for zer o-pressur e-gr adi en t flow, we 
0.01 atm for a 150 mm gauge spacing. The lower 
could have pressures one-tenth of these numbers, 
well within the state of the art. 


B. Shear-Stress Gauge 

The skin-friction coefficient will be measured with a novel 
shear-stress gauge developed at NASA Langley Research Center. 18 
The gauge consists of a 9.4 mm diameter floating disk placed in a 
housing 27 mm in diameter. The disk is connected to a four-bar 
linkage which in turn has a LVDT device with a restoring force 
motor to measure displacement. The gauge, shown in Figure 6, is 
made entirely of aluminum and has been successfully tested 
(without a heater) in the 0.3 m Transonic Cryogenic Tunnel. 
Local skin friction was measured 18 at temperatures of 100*K, 
150*K, 200*K, and 300*K, and the data overlapped reasonably well 
over the different ranges. For the present experiments, the skin 
friction will be a maximum at the forward location under 
conditions of highest dynamic pressure giving ^ 303 Pa (0.044 
psi) or a load of 0.02 N (.0049 lb.). A total of twelve gauges 
(6 each in two load ranges of 0.025 N and 0.0025 N, respectively) 
are used at the six locations shown in Figure 5. The shear- 
stress housings (without shear-stress gauges) can be used for 
boundary-layer instrumentation as described below. The housings 
are placed in the plate before the surface is finished. 


C. Hot-Wire Measurements 

The use of hot-wire anemometers in this flow environment 
presents other major problems. First, the high dynamic pressures 
may require using much smaller hot wires such as 2 micron 
diameter instead of the usual 5 micron. Second, the careful work 
of Stainback et al. 19 shows that contrary to other flows, the 
sensitivity coefficients for velocity and density are not the 
same and hence separation of the velocity, density, and 
temperature fluctuations proves to be very difficult. Stainback 
et al. 19 suggest the use of multiple-wire probes under 
simultaneous sample and hold measurements. 

Figure 7 is a sketch of a hot-wire rake that is driven by a 
lead screw within the shear-stress housing. In this way, the 
permanently-mounted housing for the shear gauges can also be used 
for the hot-wire apparatus. The lead screw will be driven by a 
stepping-motor on the back of the plate. The lead screw and the 
anti-backlash bushing will be made out of Nitronic 60 stainless 
steel that has been proven successful in the NTF environments. 
The six hot wires will be calibrated with the stepping-motor 
angular displacement. A positive backup measurement will be used 
by employing the LVDT concept parallel to the motion of lead 
screw . 
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Because of the versatility of this technique, other 
measurement instruments are used, such as the device by Ng and 
Epstein 20 which is implemented to measure freestream total 
temperature and pressure fluctuations. A total-pressure tube is 
also mounted in the shear-gauge housing. 

The use of a traverse-controlled device that would make hot- 
wire measurements over the entire plate is quite desirable when 
compared to measurements at six fixed locations. The development 
of such a device in such a severe environment is by no means 
trivial, and we continue to examine the problem. However, 
because of the extremes of temperature and dynamic pressure, the 
implementation of a boundary-layer "mouse" that can move a hot 
wire in steps of 25 pm may take longer to develop than a remove 
LDV system for NTF. 

D. Motion and Load Measurements 

An accelerometer is used to measure plate vibration, and its 
signal is correlated with the hot-wire measurements. The strain 
gauges are mounted to the back of the plate and monitor plate 
deflection under load conditions. Remote model-deformation 
instrumentation are also available. 21 These measurements are 
necessary not only for certifying the test environment but for 
safety reasons as well. 


VI. Conclusions 

It is shown that it is feasible to conduct zero-pressure- 
gradient measurements on a flat plate in the NTF under extreme 
conditions. Criteria are developed for the conduct of such 
tests . 
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